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ABSTRACT 



Context. About half of the baryons in the local Universe are invisible and - according to simulations - their dominant fraction resides 
in filaments connecting clusters of galaxies in the form of low density gas with temperatures in the range of 10 5 < T < 10 7 K. The 
existence of this warm-hot intergalactic medium was never unambiguously proven observationally in X-rays. 

Aims. We aim to probe the low gas densities expected in the large scale structure filaments by observing a filament connecting the 
massive clusters of galaxies A 222 and A 223 (z = 0.21) which has a favorable orientation approximately along our line of sight. This 
filament has been previously detected using weak lensing data and as an over-density of colour selected galaxies. 
Methods. We analyse X-ray images and spectra obtained in a deep observation (144 ks) of A 222/223 with XMM-Newton. 
Results. We present here observational evidence of the X-ray emission from the filament connecting the two clusters. We detect the 
filament in the wavelet decomposed soft band (0.5-2.0 keV) X-ray image with a 5cr significance. Following the emission down to 3cr 
significance, the observed filament is ~ 1.2 Mpc wide. The temperature of the gas associated with the filament determined from the 
spectra is kT = 0.91 ± 0.25 keV and its emission measure corresponds to a baryon density of (3.4 ± 1.3) x 10" 5 (//15 Mpc)" 1/2 cm" 3 , 
where / is the length of the filament along the line of sight. This density corresponds to a baryon over-density of pi (pc) ~ 150. The 
properties of the gas in the filament are consistent with the results of simulations for the densest and hottest parts of the warm-hot 
intergalactic medium. 
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1. Introduction 

According to the standard theory of structure formation, the spa- 
tial distribution of matter in the Universe evolved from small per- 
turbations in the primordial density field into a complex struc- 
ture of sheets and filaments with clusters of galaxies at the inter- 
sections of this filamentary structure. The filaments have been 
identified in optical surveys of galaxies (e.g. [joeveer et al.lll978t 
iBaugh et al.1l2004t iTegmark et 1012004 . but the dominant frac- 
tion of their baryons is probably in the form of a low density 
warm-hot gas emitting predominantly soft X-rays. 

The existence of this warm-hot intergalactic medium 
(WHIM) permeating the cos mic-web is predicted by numeri- 
cal hydrodynamic simulations (Cen & Ostriker 1999; Dave et al. 
1200 lb . which show that in the present epoch (z < 1-2) about 
30% to 40% of the total baryonic matter resides in the filaments 
connecting clusters of galaxies. These simulations predict that 

T < 10 7 K 
cm" 3 , which 



the WHIM spans a temperature range of 10 
and a density range of 4 x 10" 6 cm 3 to 10" 



corresponds to 15-400 times the mean baryonic density of the 
Univer se. Gas at these t emperatures and densities is difficult to 
detect (iDave et alj|2001l) . 

Several attempts have been made to detect the WHIM both 
in absorption and emission, but none of them succeeded in pro- 
viding an unambiguous detection. Marginal detection of O vra 
and Ne ix absorption features in XMM-Newton RGS spectra of 



distant quasars was reported from the vicinity of the Virgo and 
Coma clusters and was interpreted as indication of the pres- 
ence of dense and hot WHIM concentratio ns near the clus- 
ters (iFuiimoto et al.l 120041: iTakei et all 120071) A bli nd search 
for WHIM was performed by iNicastro et all (120051) . who re- 
ported the discovery of two intervening systems at z > in 
absorption along the line of sight to the blazar Mrk 421 in 
the Chandra Low Energy T ransmission Grating Spectrometer 
(LETGS) spectra. However, Kaa stra et al.l (120061) carefully re- 
analysed the LETGS spectra, confirmed the apparent strength of 
several features, but showed that they are st atistically not signif- 
icant. Furthermore, the results reported bv lNicastro et al.l (|2005) 
were neither confirmed by later Chandra LETGS observations 
nor by the XMM-Newton RG S data of Mrk 421 dKaastra et al.l 
120061; iRasmussen et al.ll2007l) . The detection of the WHIM in 
absorption with the current instruments will remain difficult, be- 
cause the expected absorption features are weak and in the large 
searched redshift space to distant blazars there will always ap- 
pear statistical fluctuations with a strength comparable to the 
strength of the absorption lines expected from the WHIM. A 
search for the WHIM from the densest and hottest parts of the 
large-scale structure filaments with the current instruments is 
more realistic in emission than in absorption. 

Searches in emission with ROSAT and XMM- 
Newton have found sev eral cand idates for WHI M fila- 

ments dKull & Bohringer! 119991; IZappacosta etafj 120021; 
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Finoguen ov et al.ll2003l :fKaastra et al. 2003). However, the asso- 
ciation of the detected soft emission with the WHIM is always 
difficult and uncertain. In some cases there is still a possible 
confusion with the Galactic foreground emission or with the 
Solar wind charge exchange emission that also produces O vn 
lines and soft emission in excess to cluster emission. 

Since the X-ray emissivity scales with the square of the gas 
density, X-ray emission studies are most sensitive to the densest 
gas concentrations of the WHIM which are, according to simu- 
lations, near clusters of galaxies. Therefore, the best approach to 
probe the physical properties of the WHIM is provided by those 
distant binary clusters between which the connecting bridge lies 
approximately along our line of sight. A filament between clus- 
ter pairs observed with such a favourable geometry has a higher 
surface brightness than filaments observed with a different ori- 
entation. By observing the variation of the surface brightness in 
X-ray images extracted in the low energy band, the presence of 
a WHIM filament can be revealed as a bridge between the clus- 
ters. This is a differential test with respect to the background and 
it is less sensitive to the systematic uncertainties. 

The most promising target for detecting the warm-hot gas in 
a filament is Abell 222/223, a close pair of massive clusters of 
galaxies at redshift z ~ 0.21 separated by ~14' in the sky (which 
corresponds to a projected distance of ~ 2.8 Mpc). Assuming 
that both clusters participate in the Hubble flow with no peculiar 
velocity, the observe d redshift difference of Az = 0.005 + 0.001 
dDietrich et al.l2 002) translat es to a physica l separ ation along the 
line of sight of 15 + 3 Mpc. lDietrich et all (120051) found indica- 
tions for a filament connecting both clusters from weak lensing, 
optical, and X-ray (ROSAT) data. They found an X-ray bridge 
with a linear extent of about 1.3 Mpc (in projection) connect- 
ing the clusters at a 5<x significance level and an over-density of 
galaxies in the filament at a 7<x significance level. In the weak 
lensing image the filament was present at a 2<x significance. 

Here we present the results of a new deep (144 ks) XMM- 
Newton observation of the bridge connecting the two clusters. 
We confirm its presence by both imaging and spectroscopy, and 
we estimate the temperature, density, entropy, and the total mass 
of the gas in the filament. Throughout the paper we use Hq = 
70 kms -1 Mpc -1 , fl« = 0.3, Qa = 0.7, which imply a linear 
scale of 206 kpc arcmkr 1 at the redshift of z — 0.21. Unless 
specified otherwise, all errors are at the 68% confidence level 
for one interesting parameter (A^ 2 = 1). 



2. Observations and data analysis 

The pair of clusters Abell 222/223 was observed with XMM- 
Newton in two very nearby pointings centred on the filament 
on June 18 and 22, 2007 (revolutions 1378 and 1380) with a 
total exposure time of 144 ks. The calibrated event files were 
produced using the 7.1.0 version of the XMM-Newton Science 
Analysis System (SAS). After rejecting the time intervals with 
high particle background, we are left with 62.6 ks of good time 
for EPIC/MOS and with 34.4 ks for EPIC/pn. 

2. 1 . Image analysis 

We coadded the final 0.5-2.0 keV and 2.0-7.5 keV images from 
all instruments and both pointings. The background was sub- 
tracted and the images were corrected for vignetting. The chip 
number 4 in EPIC/MOS 1 was excluded from the analysis be- 
cause of t he anomalously high flux in the soft band (for more 
details see ISnowden et al.1l2008l) . In order to eliminate the flux 



Table 1. The CXB components determined in a region free of 
clusters emission outside of the virial radii of the clusters. The 
unabsorbed fluxes are determined in the 0.3-10.0 keV band. / 
indicates parameters kept fixed during spectral fitting. 



Comp. 


kr/r 


Flux 




(keV / phot, ind.) 


(10~ 12 erg cirr 2 s _1 deg~ 2 ) 


LHB 


kr = 0.08' 


3.4 ± 1.1 


GH 


kr = 0.17 ±0.03 


2.9 ± 0.9 


EPL 


r = 1.41' 


22/ 


PL 


T = 0.78 ±0.16 


49 ±7 



pollution from the faint point sources due to the large wings of 
the point-spread function (PSF) of XMM-Newton, we applied 
an image restoration technique (Finoguenov et al. in prep) which 
us es a symme t ric mo del for the PSF of XMM-Newton calibrated 
by iGhizza rdi (2001). The flux from the extended wings of the 
PSF of sources detected on small scales (sum of the 8" and 
16" scale) is estimated p erforming a scale-wise wavelet analy- 
sis dVikhlinin et al.f l998). This estimated flux is used to subtract 
the PSF model prediction from the larger scales and to accord- 
ingly increase the flux on the smaller scales. The uncertainties in 
the point source subtraction are added to the total error budget in 
order to reduce the significance of the residuals. 

2.2. Spectral analysis 

We extracted the filament spectrum from a circular extraction re- 
gion with a radius of 1.6', located between the clusters, centred 
at a projected distance of 5.76' (1.19 Mpc) from both cluster 
cores. We determined the background spectrum using two dif- 
ferent sets of regions which provide a consistency check. We 
used firstly a circular region with a radius of 3.2' located at 
a projected distance larger than 1.6 Mpc from both clusters, 
which is beyond their virial radii determined using weak lensing 
(r 20 o = 1 .28 ±0. 1 1 Mpc for A 222 and r 20 o = 1.55+0.15 Mpc for 
A 223: iDietrich et al]|2005l) . Secondly, we extracted two other 
circles with 1.6' radii located on the same central chip of the 
EPIC/MOS detectors, each of these circles being at the same 
projected distance to one of the clusters as the extraction region 
centred on the filament. Using the same chip and the same pro- 
jected distance to the clusters for determining both the source 
and the background spectra ensures the best possible subtrac- 
tion of both the instrumental and X-ray background compo- 
nent. Point sources were excluded from the spectral extraction 
regions. For the identification of point sources we also used 
archival Chandra data (45.7 ks). The regions used for spectral 
extraction, for both the filament and the background are shown 
in Fig.Q] 

We subtracted the instrumental background from all spec- 
tra using closed filter observations. Since the instrumental back- 
ground level is changing from observation to observation and 
it is increasing with time over the years, we normalised the 
closed filter data for every extraction region by the 10-12 keV 
count rate for EPIC/MOS and by the 10-14 keV count rate for 
EPIC/pn. We then modeled the spectra in each background re- 
gion set with the components described below and used these 
models to subtract the cosmic X-ray background from the fila- 
ment spectrum. 

We analysed the spectra using the SPEX package 
(lKaastraetalf l996) in the 0.35 to 7 keV band for EPIC/MOS 
and 0.6-7 keV band for EPIC/pn, ignoring the 1 .3-1 .8 keV band 
because of the presence of strong instrumental lines. We fixed 
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Right ascension Energy (keV) 



Fig. 1. Wavelet decomposed 0.5-2.0 keV image of Abell 222 (to 
the South) and Abell 223 (the two X-ray peaks to the North). We 
only show sources with > 5cr detection, but for these sources 
we follow the emission down to 3<x. The filament connecting the 
two massive clusters is clearly visible in the image. The regions 
used to extract the filament spectrum and the regions used to 
determine the background parameters are indicated with red and 
yellow circles, respectively. 

the Galactic absorption in our model t o the value deduced from 
Hi radio data (N H = 1.6 x 10 20 cirr 2 . iKalberla et al.ll2005l) . We 
modelled the background with four components: with the extra- 
galactic power-law (EPL) to account for the integrated emission 
of unresolved point sources (assuming a photon index F — 1.41 
and a 0.3-10.0 keV flux of 2.2 x IP" 11 erg s" 1 cm' 2 d eg' 2 after 
extraction of point sources. lbe Luca & Molendil2 004). with two 
thermal components to account for the local hot bubble (LHB) 
emission (k^i = 0.08 keV) and for the Galactic halo (GH) emis- 
sion (kT2 ~ 0.2 keV), and with an additional power-law to ac- 
count for the contamination from the residual soft proton parti- 
cle background. The best fitted fluxes, temperatures, and photon 
indices of the background components using the 3.2' circular ex- 
traction region are shown in Table Q] 

We model the spectrum of the filament with a collision- 
ally ionized plasma model (MEKAL), the meta llicity of whic h 
is set to 0.2 Solar (proto-Solar abundances by Lodders 2003), 
the lowest value found in the outsk irts of clusters and groups 
dFuiita et alj|2007t iBuote et al.lr2004l) . However, since the aver- 
age metallicity of the WHIM is not known we report the results 
also for metallicities of 0.0, 0.1, and 0.3 Solar. The free param- 
eters in the fitted model are the temperature and the emission 
measure of the plasma. 

3. Results 

3.1. Imaging 

We detect the filament in the wavelet decomposed soft band 
(0.5-2.0 keV) X-ray image with 5<x significance. In Fig. Q] 
we show the wavelet decomposed image produced by setting 
the detection threshold to 5<x and following the emission down 
to 3cr (for details of the wavelet decomposition technique see 
Vikh linin et al.l 1998b . There is a clear bridge in the soft emission 
between the clusters which originates from an extended compo- 



Fig. 2. The spectrum of the filament between the clusters 
A 222/223 - the data points on the top were obtained by EPIC/pn 
and below by EPIC/MOS. The contributions from the X-ray 
background and from the filament to the total model are shown 
separately. 



nent. We note that no such features are detected in any XMM- 
Newton mosaics of empty fields in observations with similar but 
also much larger depths. We do not detect the bridge between 
the clusters in the 2.0-7.5 keV image. 

To verify our result, we created an image modeling the X- 
ray emission from the two clusters without a filament with two 
beta models determined by fitting the surface brightness profile 
of A 222 and A 223. We applied on this image the wavelet de- 
composition algorithm, which did not reveal a bridge like the 
one observed between A 222 and A 223. 

If the emission in the bridge between the clusters originated 
from a group of galaxies within the filament, then this group 
would be well resolved in the image. Groups of galaxies at the 
redshift of this cluster emit on spatial scales of a;0.5' as shown in 
Fig.[T]by the six extended sources detected at the Act significance 
around Abell 222/223 marked by ellipses numbered from 1 to 6. 
The filament seen in the image is about 6' wide, which at the 
redshift of the cluster corresponds to about 1.2 Mpc. 

The observed filament is by about an order of magnitude 
fainter than t h e RO SAT PSPC detected filament reported by 
iDietrich et al.l d2005h . The ROSAT detection was based on an 
image smoothed by a Gaussian with <x = 1.75' and the fila- 
ment could be the result of the combination of the emission from 
the cluster outskirts and from a few previously unresolved point 
sources between the clusters. 

The point sources (mostly AGNs) detected by XMM- 
Newton and Chandra do not show an overdensity between the 
clusters and the uniform contribution from the distant unresolved 
AGNs was subtracted from our image. Our detection limit for 
point sources is ~ 0.8 x 10" 15 erg s' 1 cm' 2 , which corresponds 
to a luminosity of 10 41 erg s" 1 for the faintest resolved galaxy at 
the cluster redshift. The number of X-ray emitting galaxies with 
this luminosity required to account for the observed emission in 
the filament would be ~ 10 times larger t han is expected from 
the distribution function o f|HasingerJ()l998), assuming a filament 
overdensity of pi (p) ~ 100. 



4 



N. Werner et al.: Detection of hot gas in the filament connecting the clusters of galaxies Abell 222 and Abell 223 



3.2. Spectroscopy 

The clusters of galaxies Abell 222 and Abell 223 are bright and 
relatively hot. The temperature of A 222, extracted from a re- 
gion with a radius of 2' is 4.43 ±0.11 keV and the temperature 
of southern core of A 223 extracted within the same radius is 
5.31 ± 0.10 keV. However, in this paper we focus on the bridge 
connecting the two clusters. 

As shown in the Fig. [2] the spectrum extracted from the fil- 
ament shows between ~ 0.5 and ~ 1 .0 ke V ~ 30% emission in 
excess to the background model (determined outside of the clus- 
ter virial radii, see Table 1). We fit this excess emission with a 
thermal model with a temperature of kT = 0.91 +0.24+0.07 keV, 
where the first error is the statistical error and the second error 
is due to the uncertainty on the background level, determined 
by varying the background components in the spectral fit within 
their lcr errors. The best fit emission measure (EM — J « H n e dV, 
where «h and « e are the proton and electron number densities 
and V is the volume of the emitting region) of the filament is: 
(1.72 + 0.50 + 0.45) x 10 65 ctrT 3 , where the first error is again 
the statistical error and the second error is due to the uncertainty 
on the background level. Assuming both errors have a Gaussian 
distribution, the emission measure of the plasma in the filament 
within our extraction region is EM = (1 .72 + 0.67) x 10 65 cm -3 
(2.6(f). Fitting individually EPIC/MOS and EPIC/pn the best fit 
temperature is 0.94 + 0.23 keV and 0.85 + 0.30 keV, respec- 
tively, and the statistical significance of the excess emission is 
3.5cr and 1.8cr. The derived emission measure depends on the 
assumed metallicity. For a higher metallicity we obtain a lower 
emission measure (for Z = 0.3 Solar, EM = 1.3 x 10 55 cm -3 ), 
while for a lower metallicity the obtained emission measure will 
be higher (for Z = 0.1 Solar, EM = 2.4 x 10 65 ctrT 3 and for 
Z = Solar, EM = 4.13 x 10 65 cm -3 ). The plasma temperature 
does not change significantly as a function of the assumed metal- 
licity. The 0.3-10.0 keV luminosity of the filament within our 
extraction region is 1.4 x 10 42 erg s _1 . In the COSMOS survey, 
a number of groups have been detected at these luminosities and 
redshifts, but their typical spatial extent (rsoo ~ 1.3', Finoguenov 
et al. in prep.) is much smaller than that of this filament. 

By determining the background properties using the two ex- 
traction regions at the central chip of EPIC/MOS and fitting only 
the EPIC/MOS data, we obtain a best-fit emission measure and 
temperature for the filament of EM = (1.7 + 0.6) x 10 65 cm 4 
and kT = 0.80+0.19 keV. These values are consistent with those 
obtained using a background region located at a larger distance 
on a different part of the detector. This minimizes chip-to-chip 
variations of the instrumental noise and allows to include some 
residual cluster emission in the form of an increased best-fit flux 
of the power-law model component. However, part of this in- 
creased power-law flux could also be due to emission from the 
filament in the background regions, resulting in the oversubtrac- 
tion of the filament. 

4. Discussion 

We detect X-ray emission from a bridge connecting the clusters 
A 222 and A 223. The temperature of the gas associated with 
this bridge is kr = 0.91 ± 0.25 keV. If this gas is the intra-cluster 
medium (ICM) at the cluster outskirts, then assuming an emit- 
ting volume with a line of sight depth of 2.5 Mpc, our best fit 
emission measure corresponds to a density of ~ 1 x 10 -4 cm -3 
and the entropy of the gas is s = kT/n 2 ^ 3 ~ 420 keV cm 2 . 
Because of the heating by accretion shocks, the cluster entropies 
rise towards the outskirts to values higher than 1000 keV cm 2 



outside O.5R200 dPratt et a l. 2006) and the entropy is expected to 
rise further out to the virial radius. The low entropy of the bridge 
implies that the emitting gas is in a filament which was not yet 
reached by the shock-heating operating on the ICM in the cluster 
outskirts. 

If the redshift difference between the clusters is only due to 
the Hubble flow, their line of sight distance difference is 15 ± 
3 Mpc. Assuming that the observed emission originates from hot 
gas within the filament connecting the clusters the mean baryon 
density of the gas is 3 .4x 1 0~ 5 (Z/ 1 5 Mpc) - 1 12 cm -3 (where I is the 
length of the filament along the line of sight), which corresponds 
to 150 times the mean baryon density of the Universe. From this 
density (assuming / = 15 Mpc) and the temperature of kr = 
0.91 keV we obtain an entropy s = kT/n 2 ^ 3 ~ 870 keV cm 2 . For 
a filament with a cylindrical shape, which is 1.2 Mpc wide as 
indicated by the images, assuming a constant density we obtain a 
baryon mass of 1.8 x 10 13 (//15 Mpc) 1 / 2 M . Assuming a baryon 
fraction of 0.16, the total mass of the filament would be « 1.1 x 
10 14 M . 

The inferred properties of the gas in the filament are remark- 
ably c onsistent with the results of the simulations by Dola g et ail 
d2006l) . They predict that within a distance of « 15 Mpc from 
a massive cluster, filaments connecting clusters have a baryon 
over-density of p/ (pc) > 100 inside a radius of «1 Mpc from 
the centre of the filament. They also found that the temperature 
of the gas along the bridges connecting some nearby clusters 
can be around 10 7 K, which is the temperature of the gas that 
we see between A 222 and A 223. The predicted X-ray s ur- 
face brightness of filaments is according to Dola g et al] d2006l) at 
most 10~ 16 erg s _1 ctrT 2 arcmin 2 . However, this predicted sur- 
face brightness was derived for zero metallicity (for a metallicity 
of Z = 0.2 Solar the predicted X-ray flux would increase due to 
line emission by a factor of ~ 1 .6) and for the case of looking 
at the filament edge on. Line emission and a favourable geome- 
try can increase the surface brightness to our observed value of 
10~ 15 erg s _1 cirT 2 arcmin -2 . 

If the observed redshift difference is entirely due to peculiar 
velocities (Az = 0.005 corresponds to Av r = 1500 km s ) and 
both clusters are actually at the same distance from us, then the 
virial radii of the clusters partially overlap and they are at the 
onset of a merger. However, if this is the case, then we do not 
expect to observe such a large amount of gas at these low tem- 
peratures between the clusters. Most of the mass from a previ- 
ously present filament would be accreted by the clusters and the 
gas in the interaction region would be compressed and shock- 
heated to temperatures significantly higher than those observed. 
The redshift difference thus cannot be entirely due to the peculiar 
velocities. 

The temperature and the average density of the observed gas 
associated with the filament indicates that we are detecting the 
hottest and densest phase of the WHIM. It is only detectable be- 
cause of the favourable geometry of the massive filament con- 
necting two large cl usters of galaxie s. However, according to 
the simulations (e.g. lDave et a l. 2001) the dominant fraction of 
the WHIM resides in a lower temperature and density phase, 
the existence of which still remains to be proven observation- 
ally. The detection of the dominant fraction of the WHIM will 
only be possible w ith dedicated future instrumentation (e.g. see 
iPaerels et al.ll2008l) . 

In order to further investigate the physical properties of the 
densest and hottest gas phase permeating the cosmic web, the 
observed sample of cluster pairs where the large scale structure 
filaments can be observed along our line of sight needs to be 
enlarged by other promising systems. 
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